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' We analytically study coarsening dynamics in a system with nonconserved scalar order parameter, 

\ when a uniform time-independent shear flow is present. We use an anisotropic version of the Ohta- 

. Jasnow-Kawasaki approximation to calculate the growth exponents in two and three dimensions: 

for d = 3 the exponents we flnd are the same as expected on the basis of simple scaling arguments, 
that is 3/2 in the flow direction and 1/2 in all the other directions, while for d = 2 we flnd an 
unusual behavior, in that the domains experience an unlimited narrowing for very large times and 
a nontrivial dynamical scaling appears. In addition, we consider the case where an oscillatory shear 
- is applied to a two-dimensional system, finding in this case a standard t^^'^ growth, modulated by 

periodic oscillations. We support our two-dimensional results by means of numerical simulations 
and we propose to test our predictions by experiments on twisted nematic liquid crystals. 
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I. INTRODUCTION 



When a statistical system in its homogeneous disordered phase is suddenly quenched below the critical temperature, 
deep into a multi-phase coexistence region, a dynamical process known as coarsening, or phase-ordering, results: 
^ ■ domains of the different ordered phases are formed and compete with each other in the attempt to break the symmetry 
and project the system on to one single equilibrium state An equivalent phenomenon occurs in the case of binary 
fluids: a system at the critical concentration tries to phase-separate after the quench, by forming domains of the two 
I different components (spinodal decomposition). An interesting problem is the analysis of the dynamical evolution 
■ O of these domains, and in particular the determination of their growth rate. In this aim a property shared by many 
^ I statistical systems, called dynamical scaling, stating that space and time scale homogeneously in the equal-time two- 
point correlation function, C{r,t) = f[r/L{t)], proves very useful. It is then natural to identify the length scale L{t) 
as the typical size of the domains during coarsening. This length scale has generally a power-law dependence on 
time, L{t) ~ t^^^, sometimes with logarithmic corrections. The determination of the exponent z for many different 
statistical systems has been the object of much effort in the past years and we can say that ordinary coarsening is 
. now quite well understood 
' A related topic, which is attracting growing attention in recent years, is the problem of phase-ordering when the 
. system is subject to an external shear. Apart from the great technological relevance of such a problem, especially in 
the case of spinodal decomposition, the basic theoretical understanding of the phenomena involved is far from being 
. well established |2|-|^]. When a shear is present, domain growth is heavily affected by the presence of the induced 
' flow and the dynamical scaling behavior is drastically different from the case of ordinary coarsening. In particular, 
"Y"; ' two main points are worthy of careful investigation: First, the growth of the domains is anisotropic and therefore the 
' dynamical evolution is described by more than one exponent. The determination of the shear exponents is, of course, 
of the uppermost importance. Secondly, it is not clear whether the shear causes an interruption of coarsening, the 
dynamical balance between growth and deformation giving rise to a stationary state (as argued in ||^), or, on the 
contrary, whether domain growth continues indefinitely. Experimental, numerical and theoretical evidence concerning 
both these points is still very tentative Q . 
^ , In the present work we perform a theoretical investigation of the coarsening dynamics in a statistical system with 
nonconserved scalar order parameter (model A, in the classification of Hohenberg and Halperin |0]), when a shear 
flow uniform in space is present. If, on one hand, such a model is unsuitable for describing spinodal decomposition in 
binary fluids, on the other hand it allows us to compute the growth exponents in any spatial dimension, in the context 
of a suitably modified version of the classic Ohta-Jasnow-Kawasaki (OJK) approximation When considering the 
relevance of nonconserved dynamics for an advance in our understanding of domain growth in the presence of a shear, 
we must take into account the fact that the only existing analytic calculations of the growth exponents for spinodal 



*E-mail:andrea@al3. ph. man. ac.uk 
tE-mail:bray@al3. ph. man. ac.uk 
■'•E-mail:rui@al3. ph. man. ac.uk 



1 



decomposition (conserved dynamics, or model B) have been performed in the limit of infinite dimension N of the 
order parameter where no saturation of coarsening is found. However, in that case the very concept of domains 
is meaningless and thus a calculation which takes into account the more physical case of a scalar order parameter 
is desirable. Besides, an understanding of the effect of shear on nonconserved coarsening is by itself an interesting 
problem, both from the theoretical and experimental point of view. Indeed, experiments have been performed in the 
past on twisted nematic liquid crystals showing that these systems are a perfect test for analytical results in 
statistical models with nonconserved order parameter. Many results, from growth laws to persistence exponents, have 
been successfully tested on twisted nematic liquid crystals ]TT[ and we therefore propose a shear experiment on such 
systems to check the results of our calculation. 

We investigate two very different cases: in the first, a shear uniform in time is applied and the behavior of the system 
is analyzed asymptotically for very large times. In the second case, we consider a shear flow which is periodically 
oscillating in time and we study the properties of the model for times much longer than the period of the oscillation. 
The primary effect of the shear flow is naturally to stretch the domains in the direction of the flow, such that they can 
be roughly represented as highly elongated ellipsoids, with the growth taking place along the main axes. Two natural 
length scales therefore arise, Ly and L±, the size of the domain along the largest and the smallest axes, respectively. 
The determination of the growth laws for these two length scales is the main objective of this work. 

In the case of a time-independent shear, our results are nontrivial and, especially in two dimensions, quite unex- 
pected. For d = 2 our calculation gives L|| ^ (Int)^/* and Lj_ ^ -y^i/^ (Ini)^^/**, where 7 is the shear rate, 
while for d = 3 we find L|| ~ ^t^^'^ and Lj_ ^ t^^^. The three-dimensional exponents are the same as one would 
expect on the basis of simple scaling arguments and are compatible with calculations for conserved dynamics in the 
large- limit Q: the growth along the flow is enhanced by a factor ^t, while the transverse growth is unaffected by 
the shear. On the other hand, the two-dimensional result comes as quite a surprise: the short size of the domains 
L± goes asymptotically to zero for very large times, while the scale area grows as in the unsheared case, L^^L± ~ t. 
As we shall show, there are topological arguments supporting this last result. As long as our approach is valid, we 
do not find any evidence of the onset of a stationary state giving rise to an interruption of the coarsening process. 
However, in two dimensions our calculation breaks down when the thickness of the domains becomes comparable 
with the interfacial width. We cannot say what happens when this stage is reached, but it is possible that some 
kind of stationary state occurs in this regime. In the case of an oscillatory shear flow in two dimensions, we find 
L|| ~ f^/^ ^/jjuj f\\{t) and L± ~ t^^^ y/^/l f±{t), where u) is the frequency of the periodic flow: both length scales 
grow like t^/^, but are modulated by oscillatory functions, /||(t) and f±{t), with the same period as the flow and with 
mutually opposite phase. In this case also, therefore, we do not flnd any stationary state. 

The structure of the paper is the following. In Section II we introduce the OJK approximation, with the appropriate 
modifications due to the presence of the shear. We end Section II by formulating some self-consistency equations for 
the matrix encoding the anisotropy of the domains (the elongation matrix). Given the technical difficulty of such 
equations, in Section III we present some simple geometric arguments useful to achieve a better understanding of the 
asymptotic behavior of the many quantities involved in the calculation. The explicit solution of the equations in two 
dimensions, together with the calculation of the growth exponents, is carried out in Section IV for a time independent 
shear rate and in Section V for an oscillatory shear rate, while in Section VI we solve the time-independent problem 
in three dimensions. In Section VII we present some numerical simulations in two dimensions, supporting our results, 
and in Section VIII we discuss a possible experimental test in the context of twisted nematic liquid crystals. Finally, 
we draw our conclusions in Section IX. A shorter account of part of this work can be found in [Of. 



II. THE OJK APPROACH 



The time evolution of a statistical system with nonconserved scalar order parameter (/>(£, t) is described by the 
time-dependent Ginzburg-Landau equation [l^ 

^ = v^0(f,.)-n0), (1) 

where V{(j)) is a double-well potential. Under the hypothesis that the thickness ^ of the interface separating different 
domains is much smaller than the size L of the domains, it is possible to write an equation for the motion of the 
interface itself, assumed to be well localized in space. This is the AUen-Cahn equation jljl, asserting that the velocity 
V of the interface is proportional to the local curvature 

v{x,t) = -V • n{x,t) , (2) 
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where n{x, t) is the unit vector normal to the interface and V • n{x, t) is the curvature. The normal vector can be 
written in general as 

^ Vmjx, t) 

Hx,t)^ ^ , (3) 

\\/m[x, t)\ 

where m{x,t) can be any field which is zero at the interface of the domain, defined by the vanishing of the order 
parameter (j)(x,t). Given that this is the only restriction on the field m{x,t), it is convenient not to use the order 
parameter itself in order to describe the motion of the interface via equation (H), but a smoother field Indeed, 
as we shall see, the principal effect of the OJK approximation is to produce a Gaussian distribution for the field 
m{x,t), which would be particularly unsuitable for the highly non-Gaussian, double-peaked distribution of the order 
parameter (j){x, t). 

From equations (||) and @ we have 

t) - - I ] _ V^m(x, t) ^ dam{x,t)dbm,{x, t) dadbm{x,t) 

By considering a frame co-moving with the interface, we can write 

dm(x,t) dm(x,t) ^ ^ ,^ , 
0= ^^+vtofym{x,t) . (5) 

If a shear flow is present, it can be taken into account by including in the total velocity Vtot of the interface, a 
contribution due to velocity field u induced by the shear 

Vtot ^ vn + u , (6) 

where vn is the curvature driven velocity, with direction orthogonal to the interface and modulus given by (jj). By 
substituting relation into equation and by noting that n ■ ~ |Vm|, we finally get the OJK equation 

a=l a,f)=l 

This is an exact relation for the field m{x,t). The OJK equation is highly nonlinear due to the dependence of the 
vector fi on the field m through eq.(||). The OJK approximation consists in replacing the factor ria^ib by its spatial 
average 

Dabit) = {naix,t)nb{x,t)) . (8) 

Note that the elongation matrix Dab must satisfy the obvious sum rule 

d 



In the isotropic case {u = 0) the elongation matrix is just Dab — Sab/d by symmetry, and the OJK equation reduces 
to a simple diffusion equation with diffusion constant equal to {d — l)/d. On the other hand, when a shear flow is 
present the matrix Dab must encode the anisotropy induced by the shear and it must therefore depend on time, as 
the average shape of the domains does. The system of equations we have to solve is therefore 

drn^ ua'-^ ^ V-Hx, t) - Y: Da^itf-^^ , (10) 

dt ^ dXa OXaOXb 

a=l a, 6=1 

n m / dam{x, t) dbmjx, t) \ 

In this paper we will consider a space-uniform shear in the y direction, with flow in the x direction. The velocity 
profile is therefore given by 
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(12) 



where 7 is the shear rate and Cx is the unit vector in the flow direction. In the present Section we will consider the 
case of a time- independent shear rate 7. A straightforward generalization of the calculation to the periodic case will 
be given in Section V. 

By going into Fourier space we can rewrite equation ( p^ as 



dmik.t) , dm{k,t) 



dt 



dky 



'kl+ J2 Dab{t)kakb m{k,t) . 



Note that a naive scaling analysis of the left-hand side of this equation would give 

Lxit)^-ft Ly{t) , 



(13) 



(14) 



where and Ly are the characteristic domain sizes in the x and y directions respectively. If we assume that the 
domain growth in the directions transverse to the flow is not modified by the shear, we obtain from (n3) the results 



(15) 



where Ly now represents any transverse direction. This is the simple scaling we mentioned in the Introduction. As 
we shall see, result (15) only holds in three dimensions, while a completely different situation occurs for d = 2. 
In order to solve equation (031) we perform the change of variables 



Qx — kx 

Qy ^ ^ ^kxt 

qa = ka , V a > 3 



(16) 



introducing the field fi^q^r) = m{k,t). The corresponding equation for fi reads 

a=3 



-qI - ily - IQxry + 

a=3 

d 

Dii{T)ql + 2Di2{T)qx{qy ~ ^Qxt) + D22{T){qy - '^qxTf + ^ Dab{T)qaqb 



(17) 



a, 6— 3 



The original OJK equation (0), with a shear flow given by (12), is invariant under any transformation which preserves 
the sign of the product xy. In order to keep this symmetry, it is necessary for the elongation matrix Dab to have the 
following block-diagonal form 

Dla{t) = D2a{t) = , Dab[t) = 2?33(i) 5ab , V a, 6 > 3 , (18) 

where, to simplify the notation, we have used D^^^it) to denote all the diagonal elements for a > 3. Equation (|lj) can 
now be integrated to give 



lJ-{q. t) = 0) exp 



-\^qaRab{T) 



qt 



(19) 



with 



i?ll(T 
i?i2(r 

i?22(T 
Rla{T 
Rab{T 



4 / dr'{[l-Z?n(T')] + 27r'i?i2(T')+7V^[l-D22(T')]} 
Jo 

4 r dr' {-D,2{t') - 7t'[1 - D22{t')]} 
Jo 

4 / dr' [1-D22{r')] 
R2a{T) =0 , V a > 3 

A5ab f dr' [1 - Dsair')] = Rsair) Sab , V a, 5 > 3 . 
Jo 



(20) 
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We can now go back to the original field m{k, t), via the relation 



m{k, t) = fi{kx, ky + jkj, k3,...,kd,t) , 



(21) 



to obtain 



m{k, t) = m{kx, ky + jk^t, fcs, . . . , kd, 0) exp 



j^kaMabit) kb 



(22) 



with 



Mn{t) = Rn{t) + 2jtRi2{t) + jh''R22{t) 

Mi2{t) = Rl2{t) + jtR22{t) 
M22{t) = R22{t) 

Mia{t) = M2a{t) = , Va>3 

Mabit) = R33{t) Sab ^ M33{t) Sab , V a, 6 > 3 



(23) 



Relation ( [2^ ) can be better understood in real space: due to the shear flow, the field m at point (x, y, . . .), at time t, is 
the propagation of the initial condition at point {x — jyt, y, . . .). Note that, if we assume a Gaussian distribution for 
m{k,0) (disordered initial condition), the field maintains a Gaussian distribution at all the times, due to the linearity 
of equation (^0|). In order to get the correlation of m{x, t) in real space we have to average over the initial conditions 

(24) 



(m(fc, 0)m{k', 0)> = V A S{k + k') 
The equal time pair-correlation function of ni is therefore 



Cm{x,x! ]t) = {m{x,t)m{x' ,t)) 



' (277)^^/2 A 

detM{t) 



exp 



].Y^ra[M-\dt) Tb 



(25) 



where Ta = Xa~ x'^^. All the information on the domain growth is contained in the correlation matrix Mabify- Indeed, 
the eigenvectors of Mab{t) give the principal elongation axes of the domains and the square roots of its eigenvalues 
give the domain sizes along these axes. 

The correlation matrix is connected to the elongation matrix by equations ( po|) and (p3|). In order to close the 
problem we have thus to write another set of equations, relating Mab{t) and Dab{t), by exploiting relation ([ill). If we 
introduce the field (fa{x, t) = dam{x, t), we can write 



Dab{t) = / VP{^) 



ifajx, t)ipb{x,t) _ 1 



II dy J VP{ip) 



e 2 



(Paix,t)(pb{x,t) , 



(26) 



and we have thus to work out the probability distribution 'DP{lp). The field tp is Gaussian and therefore we just need 
to compute its correlator. From equation (25) we have 



{^a{S,t)^b{x,t)) ^K[M-\b{t) 



'(27r)'^/2A 
dei M{t) ' 



and therefore 



'DPiifi) = — exp 



V ipa{x, t)Mab{t)(pb{x, t) 



where the constant Z normalizes the distribution. By defining 

Nab{y.t)=Mab{t)+vSab , 

and by performing the rescaling Lp — > Lp^fn^ y ^ y/n, we can write 



(27) 



(28) 



(29) 
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1 /"^ 

Dabit) = - y dy 



1 p /•D(^e-^S"'f'^'-'(^'*)^"''''(?''*)w(5,*)<^^(f,t)(^j,(f,t) 



e 2 



Let us introduce the following parameters in order to explicitly write the relation above: 

a{t) = R''hi{t)M22{t) - Ml2(i)' = Rll{t)R22{t) - Rl2{tf , 

T(t) = Mn(t)+M22(i) . (31) 

The first equation is a particular case of the more general relation det M = det R, a consequence of the fact that ( |l^ ) 
is an orthogonal transformation. We can finally write 



DlHt) = \^a{t)A-h,{tY-^ r dy 



Mll{t)+y 


[y^^ 


VT{t)y 


+ a(t)]i[M33(i)- 

Afi2(t) 


l-yj 2 




^T{t)y 


+ f7(t)]i[M33(i)- 
1 


- d-2 
l-y] 2 




VT{t)y 


+ (T(t)]^[M33(t)- 


"ryV 



2 

-^i2(t) = ^V^(i)M33(0''-' / . . . ... ... . (32) 



Relations (^0|), ( |23| ) and ( p2| ) form a closed set of equations for the correlation matrix Mabit), or, equivalently, for the 
elongation matrix Dab{t)- Before attempting to solve them, it is helpful to use physical considerations as a guide to 
the expected asymptotic form of the elongation matrix in the limit of very large times. To this aim, we will consider 
the case of a time-independent shear rate. 



III. PHYSICAL CONSIDERATIONS ON THE ELONGATION MATRIX 

When a time-independent shear flow in the x direction is present, the domains will be highly elongated along this 
direction and therefore most of the surface of the domains will tend to become parallel to the x direction for very 
large times. We thus expect the following relation to hold 

-Dii(i) = ("x^x) , t — > oo . (33) 

In the two-dimensional case, due to the sum rule (^), this relation implies 

D22it) ^ {nyUy) 1 , t^oo , d = 2, (34) 

while in dimensions d > 3 it is not a priori clear whether both Z?22 and 1^33 remain nonzero or not. The only thing 
we can write is 

D22{t) + id-2)D33{t)^l , , d>3. (35) 

With regard to the off-diagonal elements of the elongation matrix, it is not hard to convince oneself that the only 
nonzero ones can be Di2{t) = I?2i(0 = {n-xny) (see equation (p^)): indeed, due to the shear, the domains are 
elongated along two main axes which are not the {xy) axes, unless t — 00. Therefore, the sub-matrix D^^j^^\t) cannot 
be diagonal for any finite time. On the other hand, for t —i- 00 the two elongation axes become coincident with (xy) 
and thus we expect that 

£'12(0^0 , t^oo. (36) 

It is finally clear that no qualitative difference can exist between d = 3 and d > 3. Indeed, in this paper we will 
explicitly state the results only for d = 2 and d = 3. 

A useful exercise is to approximate a domain with an ellipsoid and compute the asymptotic value of Dabit) as a 
function of the main axes. We will do this explicitly in two dimensions and we will just quote the main results for 
d = 3. Let us call L|| and L± the largest and smallest axis of a two-dimensional ellipse. Besides, let 6 be the tilt 
angle, that is the angle between the x axis and the L|| axis (see Fig.l). 
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When a time-independent shear is apphed, it is natural to assume for t ^ oo: 

e^O , L\\->L^, (37) 

as an expression of the extreme elongation of the domain in the direction of the flow. We can now parametrise the 
tilted ellipse in the following way 

x{ijj) — — L|| coscj — 2^L± sincj 

y{uj) = cos a; + sinw , (38) 

with uj € [0, 27r] and where we have used the fact that 9 is very small. The average of any quantity A along the 
perimeter of the ellipse can now be calculated as, 

f!^^ duo A(uo) 

(^) = ,2. ' ' , (39) 
Jo duJ^l[u)) 



with the metric /i given by 



It is useful to compute explicitly the normalizing factor in (pQ]), i.e. the asymptotic perimeter of the ellipse, 

2 t2 /X 



^l{^JJ)^]^{L\s^T?Lo + L\cos^uf''^. (40) 



(41) 



r 1 f L L 

I duj fx{u)) ^ -L^i / duj {sin^uj^ y cos^ w)^/^ 2L|| H ^ In 

Jo 2 Jo L|| L|| 



where we have used the relation L|| 3> L±. The asymptotic perimeter divided by the total area, -L||L^, is the 
interfacial density p of the domains, which must be proportional to the energy density E of the system. In the elliptic 
approximation we therefore have 

2 

11 



i.^p^- + ^ln(^) . (42) 



It will be interesting to compare this simple result with the one obtained from the O JK calculation in the next Section. 

The vector normal to the interface can be easily found by imposing its orthogonality with the tangent vector 
{d^x,d^y). This gives 

—9LiisiniLj + L±cosu) 
^(w) 

Lii sincj + 9L i cos uj 
nyito) = J — ^ . 43 
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We can now use the relations above to compute the elongation matrix of the ellipse, Dab = (nanb)- By doing this we 
get 

Di2{t) 0^0. (44) 

Note that a priori we cannot say which one of the two pieces of Dn is going to dominate in the limit t oo. 

In dimension d = 3 it is possible to perform a similar analysis, by introducing a third axis orthogonal to the 
(xy) plane. The result is 

^2 

Dn{t) - 9^ + ^ 
^11 

Di2it)^ -6^0 . (45) 

Besides, it is possible to show that if the ratio Lj_/Lz remains constant for t oo, then both D22 and D33 are nonzero 
in this limit, and 

D22{t) = <?2(i±/i.) 

D33(t) = <?3(i±/i.) , (46) 

where the two scaling functions must satisfy the relation 

g2{x) + gsix) ^ 1 . (47) 

The results of this Section confirm our expectation on the behavior of the elongation matrix and also give us some 
hint on the relation between the elongation matrix and the domain sizes, whose determination is, of course, our final 
goal. 



IV. TIME INDEPENDENT SHEAR IN TWO DIMENSIONS 



Finding a solution of the set of equations (|2C|), (g3|) and (^ is, even in two dimensions and with a time- independent 
shear rate, not entirely straightforward. Therefore, we will first try to exploit a naive scaling analysis to find a suitable 
ansatz for the elongation matrix, and eventually we will modify our initial guess in such a way to self-consistently 
satisfy all our equations. 

First, note that in two dimensions it is relatively simple to compute the integrals in (p3). We obtain 



r(t)2-4a(t) 



where it is easy to check that sum rule (^ is satisfied. Note that, of course, eqs.(|48|) are valid also for a time-dependent 
rate ^{t) and we will therefore use them also in the next Section in the case of an oscillatory shear. 

A crucial task is now to understand which terms dominate in the limit t — + 00 in the equations above. A useful 
starting point is the correlation function in equation (|25|): if we assume that there are just two length scales, Lx{t) 
and Ly{t), a naive consequence we can draw is the following 

Mn(t) ~ L,{t)^ 

Mi2{t) ^ L,{t)Ly{t) (49) 
M22{t) ^ Ly{tf . 

Moreover, the physics of the system suggests that 

L^{t) » Ly{t) . (50) 
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Note that and Ly do not in general coincide with L\\ and Lx, as defined in the last Section. Indeed, this is the 
main difference between the naive approach and the final full solution in two dimensions. Relation ( |50| ) implies 

Mn(t) » Afi2(t) » , (51) 

and thus 

T{t) ^ Mn{t) , a{t) « Miiitf . (52) 

In order to find the asymptotic behavior of equations (|4^) we need an extra relation. From definition (^l|), it seems 
natural to assume that a{t) ^ Mii{t)M22it), and therefore, from mWj, that 



a(i) >M22(t)' . (53) 

What we are (naively) assuming is that there are no cancellations in cr(t). This assumption will fail in the final 
solution, but it will only logarithmically fail, such that relation ( |53| ) will still be true. By using relations (52) and ^2 
in (48), we finally obtain 



D,,it)^^^ (54) 

at leading order for t 00. Substituting relations (^) into (jsj) and (p5|), and by using the naive scaling relation 
Lx{t) ~ 'ytLy{t), obtained in Section II, we get 

where again we have assumed that (7{t) ^ Mii{t)M22{t). If we now use this asymptotic form of the elongation matrix 
in relations ( ^ ) and (|2^), we obtain 

Mn(0 =7'^' ^22(0 

Mi2{t) = jt R22{t) (58) 

M22{t) = R22{t) , 



and 



with 



a{t) = Rii{t)R22{t) , (59) 

i?ii(0 =472 f dt' t'^Dnit') 
Jo 

i?22(t) =4 / dt' Duit') , (60) 



always at leading order. Relations ( |58| ) are consistent with equation (pl\), and by substituting ( p8| ) into (jS^) we find 
self-consistently the asymptotic form Di2{t) ~ — l/7i. Moreover, by assuming once again that a{t) ~ Mii{t)M22{t) 
and by substituting (jSq ) into (54) we get Z?ii(t) ^ l/7i and all our assumptions seem thus to be self-consistent. 
Unfortunately, this is not the case and it is not hard to understand that something is going wrong. Indeed, if 
we now plug into equation (^) the form of (j{t) coming from equation ([59|) , rather than the naive assumption 
a{t) ^ Mii{t)M22{t), we get the following self-consistent equation for Dii{t): 

1/2 



^Rnm22{t) 1 J.dt't'^Dnjt') , 
"""^'^^ 7^t^i?22(t) J^dt'D,,it') ) ■ ^''^ 



9 



If we insert into the r.h.s. of this equation the asymptotic form of Dn found above, we find an unpleasant surprise, 
that is 



Dn{t) = ^ _ , (62) 

with a = I/V2, in contradiction with ([s^). However, the situation is far from being desperate, because if we try this 
very form of Dn in equation ( |6l] ) we fortunately find self-consistency with a = 1/2. Our initial result (|56| ) only failed 
to capture a logarithmic correction and it is possible to check that, with this new form of Dn, we recover all the 
relevant relations of this Section, namely @, (||), (H), (|55|), (||), (||), (||) and (|5l|), but not 



Summarizing, the correct final form of the elongation matrix in the two-dimensional case is therefore (always at 
leading order for t ^ 00): 



Dn{t) 



1 



D,2{t) = -^^ (63) 



D22{t) = 1-Dn{t) 
while equation (^6|) is not correct. From ( |60| ) we have 

7 



i?n(t) - (64) 



R22{t) 



4 Vln 7^ 

7 



whereas, from (pq), the correlation matrix is 



Mii(t) = 47<2 0n7t 

Mi2{t) = At ^/h^t (65) 

M22{t) 



4 Vln 7t 



7 

cr(t) = 4^2 . 

It is possible to see now that the critical assumption which went wrong in our initial analysis was a{t) ~ Mii{t)M22{t)- 
Indeed, from equations (65) we see that CF{t) is smaller than this, because there are some non-trivial cancellations 
in the determinant of Mah- For this same reason, one should not be misled by the fact that apparently in ( p5[ ) the 
determinant of Mab is null: we did not write the sub-leading contributions to the correlation matrix which make 
(j{t) r^t^^t^ Int. 

In order to obtain the domain size along the principal elongation axes, L\\{t) and L±(t), we have to find the 
eigenvalues \i{t) and \2{t) of Mab{t)- This is easily done by recaUing that the characteristic polynomial is just 
— rA -I- cr, where t and a are the trace and the determinant of Mab{t), respectively (cfr. eq.(|3l|)). The final result 
for the two-dimensional case is: 



til (!) = -/W) = VW) = 2^1 (In jt)'/' 




Note how striking the effect of the shear is in two dimensions: the size of the domains along the minor axis shrinks 
to zero, even though very slowly, for t — > 00. The asymptotic effect of this unlimited narrowing of the domains for 
very large times is still unclear to us. However, we do expect our approach to break down when L±(t) becomes of the 
same order as the interface thickness ^, when equation (||) ceases to be valid. This happens after a very large time, 
of the order exp(l/7^^^). What we can say is that, if a steady state exists, it can be reached only when the thickness 
of the domains becomes comparable with the interface width. 

An important feature of the solution we have found is the failure of standard {x, y) scaling. In order to appreciate 
this fact, we have to remember that, even though L|| and L± are the natural domain sizes along the eigen-axes of the 
correlation matrix, other length scales can be defined, as shown in Fig. 2: 
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FIG. 2. The length scales L^, Ly, and H, 



First of all, we have and Ly-. from the correlation function (^5|), it follows that 

1 



Ly{t) = ^^^]_^^ , (67) 



and from ( |65| ) we get 



(ln7t)i/ 



Secondly, we can define and iJj,, as the maximum extension of the domain in the x and y directions, that is 

Hx = L|| cos 6 

Hy = L\\ sine, (69) 
where 9 is the usual tilt angle (see Fig.l), which can be easily computed from the eigenvectors of M. These are 



, e^=(-\,l] , (70) 



and therefore 



In this way we have 



(71) 



7< / \ 7^ 

It 

(72) 

In the absence of shear all these length scales would coincide, that is we would have L\\ = — and L±^ ~ Ly ~ 
Hy. With the shear this is no longer true, simply because A/12 7^ 0. Still, we would expect these length scales to differ 
only by some constant factors, such that they would all be of the same order asymptotically in time. If this situation 
held, we would have a standard {x, y) scaling, even though with anisotropic domains. However, in two dimensions the 
situation is very different, because the length scales above differ by logarithmic corrections. More precisely, we have 



1 



L^[t) ^ Ly{t) ^ ^_ Hy{t) . (73) 



The fact that L\\ 7^ L^, and therefore the emergence of a nonstandard dynamical scaling, is closely related to the 
vanishing of the determinant of M at the leading order, and its consequence is that {x, y) are not the correct scaling 
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axes. We shall see that this fact does not happen in three dimensions. In order to obtain the right scaling, we have to 
refer to the eigenvectors of the correlation matrix M, from which we can finally write the scaling form of the two-point 
correlation function in two dimensions 



hit)' L±{t) 



with 



s = X + y/jt 
u — y — x/^t 



(74) 



(75) 



In the expression above, / is a scaling function, while s and u are coordinates along the main scaling axes of the 
domains. Note that by x and y we actually mean rx and Vy. 
Furthermore, note that the elongation matrix can be written as 



Dii{t) 



Li 
ill 



to be compared with the result for Dab obtained with the elliptic approximation (cq.(}44D). 

An interesting quantity which can be easily computed is the interfacial density p(t), defined as 

p{t) = ( S{m{x, t)) I Vm(f , t)\)= I VP{m, ^) 5{m{x, t)) \0{x, t)\ , 



(76) 



(77) 



where, as in Section II, we have put (fa{x, t) = damlx^t). The calculation is easy to do because the Gaussian fields 
m and ip are uncorrelated. From relations (Eq) and (pq) we have 



Pit) 



a{t) 



1/2 



(27r)3/2 



e 



\mt)\ ■ 



By using the following formula 



f°° _dV_ ( f 
Jo \^ 



-V y 



r ^ i^-' - 1) r(-i) 



dy 

1,3/2 



{^e-'^'y - 1) 



we can perform the Gaussian integral over ip in ( [78| ) and, by proceeding as at the end of Section II, we get 



Pit) 



V2 



dy 

r(_i)(2vr)3/2 7o 2/3/2 



(1 + yr/a + y'^/a)^/'^ 



In 



(78) 



(79) 



(80) 



where we have used the asymptotic expressions (^ for T(t) and a(t), together with relations (|6|). Remarkably, this 
formula for the interfacial density has the same asymptotic form as the one we have obtained in the context of the 
elliptic description of domains (see eq.(p^)). Besides, we note an important point: p{t) is proportional to the energy 
density of the system and therefore, given that L±{t) decreases with time (equation (|66|)), equation ( |80| ) means that 
the energy in the two-dimensional case increases with time 

Eit)^-^^^ilnjt)'/\ (81) 

where we have subtracted the trivial ground-state contribution. This may seem a surprising result, but we have to 
remember that due to the shear the system is not isolated, and therefore the dynamics is not a simple gradient descent 
(in other words, no Lyapunov functional exists). A simple example can make this point clearer. Imagine we prepare 
a two-dimensional system between two boundaries in a striped configuration (see Fig. 2), with the stripes orthogonal 
to the boundaries (assume fixed boundary conditions according to the stripes). This configuration is stable at T = 0. 
If we now shear this system, by moving the boundaries in opposite directions, the stripes will be stretched and the 
interfacial length per unit area will increase (see Fig. 3). Thus, in this simple case, the energy of the system increases 
under the application of a shear. This example shows that there is no general reason why the energy of a sheared 
system cannot increase with time. Of course, it would be important to test equation (|8^), together with all our 
predictions, in a numerical simulation or even better in a real experiment (see Section VIII). 
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FIG. 3. Increasing the energy of a system by shearing it. 



The OJK theory also gives an explicit expression for the scaling form of the correlation function 
follows from equation (p5|) and from the scaling relations above 



CojK{x,y;t) = - sin 



(m(l)TO(2)) 

(m(l)2)l/2(TO(2)2)l/2 



2 . 

= — sm 

TT 



cxp 




I, which simply 



(82) 



It has been noted in p5[ that in an unsheared, but anisotropic system the OJK form of the correlation function fits 
very well the numerical data. Note, however, that in the present case, unlike in the scaling laws along the two 
main directions are radically different due to the shear, and therefore it is not a priori clear to what extent ( p2| ) is a 
good approximation to the scaling function / in ([z^). On the other hand, we believe that the scaling form we find 
in ([tJ) has a general validity. Finally, let us note the elliptic symmetry of the OJK correlation function, which could 
explain the partially correct results we obtained by approximating the domains with ellipses. The same will be true 
in three dimensions. 

An important property of the result we have found is that the scale area of a domain A{t) satisfies the following 
relation: 



A{t) = LnL^ = 2t , 



(83) 



as in the case where no shear is present. As we are going to explain, there are topological reasons why in two 
dimensions relation ( ^3|) must be satisfied either with or without shear. Equation (|8^ ) is thus a necessary condition 
fulfilled by our result, which, by itself, clearly shows that the transverse growth must be depressed if the longitudinal 
one is enhanced. 

Let us consider an isolated domain in two dimensions in the absence of shear. The rate of variation of the area 
enclosed in the loop is 

dA{t) 



dt 



— <f> dl V = — (t dl V 



(84) 



where v is the velocity of the interface and V • n is the local curvature (see eq.(||)). By virtue of the Gauss-Bonnet 
theorem, the right-hand side of equation ( p^ is in two dimensions a topological invariant, and therefore independent 
of the shape of the domain. 

When a shear is present, we have to add to the velocity due to the curvature the flow velocity u in the direction 
orthogonal to the interface. The right-hand side of equation (p3) is thus corrected by the following term: 



dl ft ■ 







(85) 



the final equality holding for any divergence- free shear flow. Equation (p3D, therefore, holds in two dimensions 
irrespective of the presence of the shear. It is interesting that the OJK approximation, in the self-consistent anisotropic 
version we have presented here, is able to capture this essential topological feature of phase ordering in two dimensions. 
Note also that the constant 2 in relation ( |83| ) is exactly the same as one would obtain from the domain size in the 
absence of shear. We will find the same constant in the case of an oscillatory shear, as a further confirmation of the 
validity of our method. 
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V. OSCILLATORY SHEAR IN TWO DIMENSIONS 



The rather surprising rosuhs we have obtained in two dimensions could raise the question whether the OJK method, 
in the modified form we are using here, is actually suitable for studying the physics of a sheared system. Indeed, 
the skeptical reader may very well think that the shrinking of the transverse domains size, with the consequent 
increase in the total energy of the system, could be an artifact of the technique, rather than a genuine property of 
the model. On the other hand, as we have seen at the end of the last Section, our two-dimensional result satisfies 
the highly nontrivial topological relation on the growth of the scale area, eq. (p3[), supporting the validity of our 
findings. Therefore, to check how robust our method is, we test its compatibility with the two-dimensional topological 
constraint in a completely different situation. To this end we study in this Section the effect of an oscillatory shear 
on phase ordering in two dimensions. 

It must be said that the case of oscillatory shear is interesting in itself. Indeed, a realistic experimental situation is 
very unlikely to involve an indefinite time- independent shear. More reasonably, a shear flow periodically depending on 
time, typically with some random modulation, is what we expect. Of course, real experiments with time independent 
shear can be performed (and we are proposing one in Section VIII) : what we are saying is that a generalization of our 
calculation to a time-dependent oscillatory shear can shed some light on a more natural experimental setup. 

We consider a sheared system with a velocity profile given by 

u = jyG{t)e^, (86) 

where the only assumption we make on the shear function G{t) is that it is a periodic function with fundamental 
frequency uj and zero time average. One of the interesting aspects of the following calculation is that the results are to 
a great extent independent of the explicit form of G{t). The derivation of the OJK equation is completely analogous 
to the one in Section II, and it follows simply from the obvious substitution 

dmjk, t) dm.{k, t) 

In order to solve the equations we have therefore to perform the change of variables (compare with (|1 

1v = + —9{t)kx , (88) 

with 

g{t) = u f dt' G(t') . (89) 
Jo 

Of course, for G{t) = 1 we reproduce the time-independent shear case. All the equations of Section II can now be 
generalized to the oscillatory shear case by means of the following trivial substitution 

7t ^ ^g{t) . (90) 

UJ 

A critical issue to understand concerns the regime of the parameters, in particular time, that we have to consider. 
First of all, we cannot afford to have too high a frequency, otherwise there would be a delay in the response of the 
system to the shear. This means we must take the shearing frequency, w, much smaller than the shear rate, 7. On 
the other hand, we need to observe the system on time scales much larger than a period. Therefore, we will consider 
the following regime 

- < - < t , (91) 

7 U! 



which implies 

_ 7 



> 1 . (92) 



Note that, in this way, we cannot recover from our final results the w — *■ case, nor can we extrapolate to the lo - 
limit. On the other hand, the large parameter a will be useful for extracting the leading terms from our results. 
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Before going further, let us explain our general strategy. Due to the periodic shear, all our quantities will exhibit 
oscillations: some of them, like Dn, which is positive definite, will oscillate around a nonzero value, while others, like 
Di2, will oscillate around zero, due to the oscillation in the orientation of the domains. Given that all these quantities 
enter the time integrals in equations (|20|), a natural approach, for times much longer than the period, is to exploit 
their time-average: if B(t) is an oscillatory quantity we write, to leading order for t — > oo, 



t 

dt' B{t') t'" - B , (93) 







with 



I dt' Bit') . (94) 





In this way from equations ( ^ ) we get 

i?ii(t) =4t{l-ri7+2a^^+a2 '^^I'j 

i?i2(i) =4t{-^-a (95) 

Note the striking difference from the time-independent shear case: due to the oscillations the whole matrix Rab is now 
of order t, as it would be in the absence of shear. As mentioned above, we expect D12 to oscillate around zero with 
the same period as g. Thus, in the equations above we can disregard terms like £'12 and gDn, whose time average is 
zero. As a consequence, we have 

Ruit) = , (96) 
that is, the isotropy is restored at the level of the matrix Rat- On the other hand, we have to keep mixed terms like 



gDi2, because their time average will be nonzero. Using equations (^3|) we can now write 

Mii{t) ^ U ^^l -D^ + 2a'^D^ + g^Dn + a^g^{t) D^} 

Mi2it)^Ua g{t)D^ (97) 
M22{t) ^UThl . 

First of all note that, apart from the oscillation induced by the explicit presence of g{t)^ the correlation matrix Mab 
is of order strongly suggesting that we will end up with a t^/"^ growth. On the other hand, M12 7^ 0, meaning that 
the system is still anisotropic, even though the anisotropy has zero time average. 

From relation (p2[) and from equations (p7|), we have that Mn 2> M12 ^ M22, and therefore the self-consistent 



equations (48), become 



I,^^^^^M22[t) + ^) 



M^2{t) 
Mnit) ■ 



As usual, we need a starting point to break into these equations and some physical considerations may help here. 
First, note that naively M12 ~ L^Ly '-^ t, from the topological relation (|83|). The second of equations ( p7| ) then 
suggests that Dn ^ \/a. Secondly, from the form of the velocity profile, we have another naive relation, that is 
X ~ ayg{t). Thus, we expect that —D12 ^ 6 ^ y/x ^ \/a. We therefore make the following ansatz 

Dii{t) - -f{t) 
a 

-Di2{t) = -h{t) , (99) 
a 

with f{t) > 0, while we expect h to oscillate around zero. Both / and h must now be determined self-consistently. 
Inserting this ansatz into (p^) and considering only the leading terms in a, we have 
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with the two constants u and r given by 



r 



Mii(i) = Ua[u + rg^{t)] 
Mi2{t)=Atg{t)r 

M22{t)=4t- (100) 
a 

a{t) = 16t^ru, (101) 
7 u = '^. (102) 



By inserting this form of Mab into equations ([O^) , we find that the powers of a balance and we obtain two equations 
for the functions / and h 



u + rg''(t) 

m^n{r/u), (104) 



Averaging the first equation, we get 



with 

On the other hand, by multiplying the same equation by g^{t) and averaging again, we have 

y/Tu = 1 - n{r/u) , (106) 

and therefore 

v^ = r2(r/u) = 1/2 . (107) 

In order to compute the domain sizes we can use the same formulae as in Section IV, because we still have ~ 
Mfi > a. We obtain 

'7 



L||(i) = V^=2ti/2 ^1 ^u + rg^t) (108) 
and, happily, we find for the scale area 

A{t)=L\\{t)L^{t)^At^fF[1^2t , (110) 

independent of the explicit form of the shear function G{t). Note also that the factor 2 in this formula is exactly the 
same as in the time-independent shear case and in the unsheared case. This is an important result, supporting the 
validity of our method for the study of the effect of shear in this type of system. 

As expected, apart from the oscillations, the growth follows a i^/^ law. The interesting thing is that both Ly and 
oscillate in time, but, as expected, with an opposite phase: when g(t) has its maximum (i.e. at the maximum 
shear displacement), Ly is maximum and of course is minimum, because this is the point of maximum elongation 
of the domains. On the other hand, for g(t) — (i.e. zero shear displacement) L|| is minimum and L± maximum, 
but always with Ly ^ L_l. We want to stress that this oscillatory dynamics is only deceptively simple. To better 
appreciate this fact we have to compute Lx and Ly (see Fig. 2 and eq.(|67|)). These quantities read 
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=2^1/2 ,p (111) 
y UJ 

Ly{t)=2t^'\l^ ^ ^ . (112) 

First of all, note that Lx, unlike Ly^ does not oscillate in time, and this had to be expected from its very definition (see 
Fig. 2). Secondly, note that for g{t) — 0, we have Lx = i|| ^ Ly = L±: at the points of zero shear displacement the 
domains are very flat and large. Besides, we can compute the tilt angle 9 from the eigenvectors of M, thus obtaining 

tan0 = i^^^%^. (113) 

We can see that is zero at the zero displacement point {g{i) = 0), and increases with increasing displacement, up 
to a maximum, whose value decreases with increasing shear rate 7. This fact may seem counterintuitive, especially 
because in the case of a time-independent shear rate we have seen that the tilt angle was decreasing with time, while 
here it is increasing. However, there is no contradiction: the key point is that at g{t) — the domains are already very 
elongated, that is Lx ^ Ly, as an effect of the shear experienced in the former periods. We can better understand 
what happens by using the simple case of a linearly sheared ellipse (no growth), with initial axes Lx and Ly, and 
Lx ^ Ly. The ellipse is described by the parametric equation 

X — Lx cos (j) + jyt 

y = LyS'mcf) , (114) 

with 4> G [0, 27r]. We can estimate the tilt angle by computing the ratio y/x at the point where the x displacement is 
maximum. This gives 

This function has a maximum at 

imax = - 7^ , (116) 

7 Ly 

and decreases asymptotically like for t 00. In the case of time- independent shear rate, the initial condition 
t = has Lx/Ly ~ 1, and therefore the maximum of 9 is quickly reached at tmax ~ I/7, which is much smaller than 
the times we consider, t 3> I/7. For this reason, in the time region of interest the tilt angle monotonically decreases. 
In the oscillatory shear case, on the other hand, at the zero displacement point, g = 0, we have Lx/Ly ~ 7/0-1 and 
thus tmax ^ I/lu: thc tilt angle therefore increases during the period of the oscillations, and this explains the apparent 
contradiction between the two cases. 



From the tilt angle (113) we can compute the additional length scales Hx and Hy by using definition (|69[). We have 



y UJ 

- 9/1/2 rgjt) 

7 ^u + rg^{t) 



(117) 

i?,(t)=2tV2J^_^-v__ (,,8) 



Note that Hy{t) is the only length scale to vanish at the zero displacement point. After the discussion above, the 
reason for this should now be clear. 

In order to compute r and u, we need to know the explicit form of g{t), and therefore of G{t). However, these 
are just numerical constants and the time evolution of the domain sizes is not afi^ected by them. For the particularly 
simple case where 

G{t)^cosujt, (119) 

the constants are 
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VI. TIME INDEPENDENT SHEAR IN THREE DIMENSIONS 



In dimension larger than two it becomes very difficult to explicitly compute the integrals in equations (^2|) . Notwith- 
standing this, if we formulate a suitable ansatz for the elongation matrix Dab{t), we can then find Mab{t) from equations 
( pO| ) and (p3|), and finally obtain a self-consistent relation for Dab{t) by an asymptotic evaluation for large times of 
the integrals in (p2|). In the present Section we will carry out this program for a time-independent shear rate. 

First of all, we note that many of the terms in equations ( pO| ) can be estimated by means of the following reasonable 
ansatz: 



Du{t) ^ -1 . (121) 

Both these relations are also obtained in any dimension by the calculation of Section III. By inspection of equations 
( pO| ) it is now clear that the key quantity needed to evaluate Rab{t), and thus Mab{t), is [1 — £'22 (^)]- We could be 
tempted to try an ansatz similar to the case d = 2, by taking [1 — D22{t)] ~ (in jt)"''^ / f^^ 0, for t 00. However, 
a careful analysis of the equations shows that this ansatz is not consistent. Therefore, the most natural thing to do 
is to assume that both D22{t) and I?33(i) remain nonzero for t — s- 00, that is (according to the usual sum rule) 

D22it)^l-K (122) 
^33(0 K , 

and to fix self-consistently the value of the constant K. From equations (pO|), (p3|), ( ^ and ( |122| ) we have 

Mi2{t) = 2K-ft^ 

M22it)=4.Kt (123) 
Masit) = 4(1 - K) t 

o{t) = ^^KWt' , 

at leading order for t 00. Note that the explicit forms of I?ii(<) and Di2{t) do not enter in Mab{t)- Using relations 
(123) it is now possible to evaluate the asymptotic value of the integrals in equations (|3^) and get an equation for the 
constant K. In three dimensions equations ( p^ ) read 



2"^ " ''"J„ ^ {y3 + Mn{t) 2/2 + [Afii(i)M33(t) + <j{t)] y + a{t)M33it)V^^ 

n (f^-}- I m r ^ y^ + M^,{t)y + a{t) 

^33W 2^^WiW33Wy^ {y3 + Mn(t)y2 + [Mn(t)M33(0+^(i)]2/ + ^(i)M33(t)P/2 



where we have used th e rela tion Mii(i) 3> M22(t) ^ M^^lt), according to (12S). By performing the rescaling y ^ ty 
and by using relations (122) in the two integrals above, it is possible to see that in the limit i — > 00 we can disregard 
the terms y^ at the numerator and at the denominator. In this way we obtain 



D22{t) = ^ dy 



2 



with 



^33(i) r dy , .3,, , (124) 

2 Jq [y^ + y + ayi^ 



4K{1 - K) 



(4-3X)2 
K 
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The fact that there is no time-dependence left in the right-hand sides of equations (124) shows that ansatz (121) and 
(122) give rise to a self-consistent solution for the three-dimensional case. Moreover, it is straightforward to check 
that sum rule (||) is satisfied. The integrals in (124) can now be easily performed and, by using relations (122), after 
some algebra we find 



A similar treatment of the integrals in 



X = 1/5 . 
|) for -Dii(i) and Di2{t) shows that. 



(126) 



72*2 
it 



consistent with ansatz (121). Let us note that relations (121), (122) and (123) are self-consistent in any dimension 
d > 3, as can be easily verified by using these relations in equations ( |3^ ) and rescaling y ^ ty in the integrals. Our 
final result will therefore be qualitatively the same for any dimension d > 3 (for d > 3 only numerical factors, such as 
the values of K and the amplitudes in equations (127) below, are changed). 

We can now compute the eigenvalues of the correlation matrix Mab(t), in order to find the sizes of the domains 
along the principal elongation axes. From (|123|) and (126), we have 



whose corresponding eigenvectors are 



^ii^(^'2:^'« 



L\\{t) 
L±{t) 

e± = 



2 


-ft" 




1 


^1/2 




4 


il/2 , 


7t 



(127) 



3 

2^' 



1,0 



(0,0,1) 



(128) 



where we recall that L|| and L_l are the larger and smaller orthogonal axes of the domain in the {xy) plane, whereas 
Lz is the axis of the domain in z direction (or any direction orthogonal to the {xy) plane, if d > 3). The domain 
growth in dimension d > 3 is therefore the one wc would expect on the basis of the simple scaling arguments given 
in Section II (see equations (p^): the growth exponent along the flow direction is augmented by one, whereas the 
others are left unchanged. Unlike the two-dimensional case, there are no topological restrictions on the product of 
the domain sizes, because the integral over the domain surface of the local curvature is not, in d 7^ 2, a topological 
invariant. 

As already anticipated, for d = 3 standard scaling holds. Indeed, one can immediately check that 



L±{t) 



Ly{t) 



Hy{t) 



(129) 



There is therefore no real difference between growth along the principal axes of the domains and growth in the {xyz) 
directions, and the correlation function displays the simple asymptotic scaling form 



C{x,y,z\t) = f 



xyz 



According to the OJK theory we have, for t ~^ 00 

2 . f 



CojK{x,y,z;t) 



exp 



1 ({x + 32//27t)2 ^ {y - ix/2^tf z^ 

LI 



(130) 



(131) 



In the scaling limit, where x,y, z,t 00 with a;/L||, y/Lj_ and z/L^ fixed, the term 3y/2jt can be dropped, but the 
term Sx/2^t cannot, and the OJK scaling function has ellipsoidal symmetry as expected. 
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As in two dimensions, we can compute the interfacial density by applying equation (fTn). The final result is 



Pit) -JT + IJ^ (132) 



which shows that the energy density in the three dimensional case decreases in the standard way, 

E{t) ~ _ (133) 

VII. NUMERICAL SIMULATIONS IN TWO DIMENSIONS 

In the present Section we will present some numerical simulations for a two-dimensional system subject to a time- 
independent uniform shear. We have considered a system of Ising spins on a lattice, governed by zero-temperature 
Monte Carlo dynamics. As in the rest of this paper, the shear flow is applied in the x direction, according to the 
profile given by equation (p^. From a practical point of view, we have sheared the system by shifting each row of 
spins by an amount proportional to the y coordinate and to the time t (measured in Monte Carlo steps) 

Aa:(y,£) =2/n,(f) , (134) 

where ns{t) is the number of discrete shear steps up to time t. Of course, the discrete nature of the system is reflected 
in the discrete nature of the shearing process. To simulate a shear rate 7 (defined by Aa;continuous — lyt)) we require 
ns{t) — Int {"it): After each I/7 Monte-Carlo steps a discrete shear process, where each row moves one lattice spacing 
relative to the row below it, is applied. In the large-time limit, where i ^ I/7, the system's behavior should not be 
very different from that of a continuously sheared system. 

We have to be careful in choosing the boundary conditions for a sheared system, because normal periodic boundary 
conditions would clearly be wrong. The idea is to replicate the original system infinitely many times on the (x, y) 
plane and to shear each sub-system with respect to the others. In other words, if {x, y) are the coordinates on the 
infinite plane, and are the coordinates on our numerical system, we have 

i = ModN^{x -\- yns{t)) 

j = Modw^(y), (135) 

where Nx and Ny are the sizes of the numerical system in the x and y directions, and the function ModAr(z) is just 
the value of z modulo N . Clearly, for i = Q equations ( |135D reduce to standard periodic boundary conditions. 

One of the main difficulties in simulating a system subject to a shear is that the domains grow very quickly in the 
direction of the flow, soon reaching a size comparable with the size of the system. On the other hand, as we have 
seen, we expect the growth to be highly depressed in the transverse direction. Thus, the most reasonable thing to do 
is to take ^ Ny, in order to reduce as much as possible finite size effects. In all our simulation we have taken 
Nx — 20000 and Ny — 100. As we shall see, even for our longest times, the domains are much smaller than the size of 
the system in both directions. A possible proposal in order to reduce the finite size effects due to the shear-induced 
elongation of the domains is to work at very low 7. However, all our results hold in the limit Lx S> Ly-. if we decrease 
the shear we will have to wait for a longer time to enter the asymptotic regime of interest, and thus we will still have 
the problem of long domains compared to the system size. There is, therefore, no easy way out of this situation and 
we had to tune our parameters to take this problem into consideration. For this reason we run our simulations for 
only one value of the shear rate, namely 7 = 1/4: in order to study the dependence of all the observables on the shear 
rate we would have to consider values of 7 far from the suitable numerical domain. 

The first thing we want to check is the behavior of the length scales Lx{t) and Lyit). As we have seen, x and y 
are not the correct scaling axes, but we want to test our prediction for Lx and Ly against the naive expectation of 
equation (^5|). Indeed, it must be remembered that this naive scaling is also the one found in the case of conserved 
dynamics in the limit of infinite dimension of the field We recall our analytic prediction 



Lx{t) 



(ln7t)i/4 
1 
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Note that Ly is, at the leading order, equal to and therefore we can limit ourselves to measure the former. This is 
important, because a numerical measure of L_l would be very difficult: the domain size in the perpendicular direction 
is very small and for long times this direction passes through very few lattice sites, such that practically there are no 
points where the correlation function is different from zero. This problem does not exists for the correlation in the 
X, y and parallel direction. In order to extract the the domain scale at a given time we have performed a fit of the 
correlation function to the OJK form and have located the point where the fit is equal to 1/e. We have checked that 
the behavior of the domain size with time is almost entirely insensitive to the particular fit we use. Numerically, we 



do not expect to be able to detect the logarithmic corrections in (136), so our goal is to check the leading behavior 
Lx ^ t and Ly ~ 0(1). Our results are shown in Fig. 4. 
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FIG. 4. The domain sizes in the x and y direction as a function of time. The full line is a power fit giving Lx ^ t ' . The 
dashed line is L"^'"" ~ t^^'^ , and the dotted line is LJ)'"^'^ ~ t^^"^ , for a comparison. In both cases, the data are averaged over 5 
samples. 



As we can see, is definitely not growing like t^^'^ . A power fit gives 

LS) - t' "' ■ (137) 

Furthermore, Ly is, on this scale, compatible with a constant, and is certainly not growing like Both Lx{t) and 
Ly(t) have the expected behavior, apart from the logarithmic corrections, and the naive exponents 3/2 and 1/2 are 
clearly not correct. 

The value of Ly is very small and in order to have a better idea of the fast decay of the correlation in the y direction, 
we plot in Fig. 5 the correlation function and the OJK fit for a given fixed value of the time. Note that actually the 
correlation vanishes on average after six lattice spacings. 

The next important quantity we want to measure is the energy. From relation ( |8l| ) we can see that E{t) is a direct 
measure of L±{t). Note that the relation between E and Lj_ is, at leading order, completely independent on the OJK 
approximation we are using: indeed, the simple assumption L|| 3> L±_ is sufficient to conclude that, at leading order, 
E ^ l/L±. However, we stress that the condition L|| ^ L± is only satisfied for large times (see Fig. 4). In Fig. 6 we 
plot the energy as a function of time, both for the sheared and the unsheared case. We see that, after an initial drop 
in the time regime where we do not expect relation (^ij) to hold, the energy becomes compatible with a constant on 
this scale. The difference with the unsheared case is striking. In the inset of this figure we show a magnification of 
the last part of the curve for the sheared case: it is encouraging to see that, despite the significantly large error bars, 
an increase in the energy for very large times is clearly visible, compatible with our analytic prediction 

Eit)^ ^i\n-ft)'/\ (138) 

However, we stress that longer simulational times and larger system sizes are necessary to test this prediction (in 
particular, the power of the logarithm) more carefully. 
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FIG. 6. The energy as a function of time in the sheared case (7 = 1/4), averaged over 41 samples, and the unsheared case 
(7 = 0) averaged over 5 samples. The full line is a power fit for the unsheared case, giving Eun ~ t'"'^^ . The horizontal broken 
line is a guide to the eye. Inset: energy in the sheared case as a function of time (magnification). 



The last quantity we measure is whose form ( pq ) differs from that of Lx{t) only by a logarithmic correction. 

In Fig. 7 we plot L|| as a function of the time. Even if slightly faster, the growth of the domains in the parallel direction 
is compatible with t. Indeed, a power fit gives 

L\\{t) ^ t^-^^ . (139) 

Not surprisingly, at a simulational level we are unable to detect any significant difference between the growth of 
and L||. 
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FIG. 7. Domain size as a function of time, averaged over 10 samples. The full line is a power fit, giving ~ t^'^*. 

Summarizing, we can say that, up to the simulational times we were able to reach, numerical data are largely 
compatible with our theoretical results. In particular, the nontrivial leading behavior of Lx{t), Ly{t) and L\\(t) are 
correctly reproduced, while the naive expectations for the domain growth is sharply ruled out by the simulations. 
Note that, of course, a longer simulational time would be desirable, especially to check whether the curve of the 
energy develops a well defined minimum and eventually starts increasing as (Int)^/'*. Unfortunately, as we have seen, 
the high values of and L|| make this impossible with the system sizes we were able to reach, otherwise finite size 
effects would heavily come into play. Also for this reason, in the next Section we propose a real experimental test of 
our analytical results. 



VIII. AN EXPERIMENTAL TEST OF THE TWO-DIMENSIONAL RESULTS 



Theoretical and numerical results on nonconserved two-dimensional coarsening dynamics can be experimentally 
tested by means of thin films of uniaxial twisted nematic liquid crystal (TNLC) subjected to rapid thermal quenches. 
Since the classic experiments of Orihara, Ishibashi and Nagaya , showing that a dynamical scaling compatible with 
the law L(t) ~ t^^'^ actually takes place in this system many other workers have successfully tested numerical 
and theoretical results on nonconserved coarsening in TNLC |11|. In particular, let us note that this kind of system 
seems to be particularly suitable for testing our analytic calculation: indeed, it has been shown in [Tc| ] that the scaling 
function describing the two-dimensional coarsening dynamics in TNLC is very well approximated by the analytic 
expression given by the OJK theory Moreover, it has been explicitly checked that the AUen-Cahn equation 
(^, describing the motion of an interface due to its curvature, holds to a very good degree of accuracy for TNLC. 
Our aim is to describe in this Section the basic experimental setup for TNLC and to propose a shear experiment on 
such systems, in order to test our nonstandard two-dimensional results in the case of simple, time-independent shear. 

A typical TNLC cell is obtained by confining the sample of nematic liquid crystal between two glass plates, previously 
prepared by rubbing them in two mutually perpendicular directions. In this way the orientations of the crystal 
molecules belonging to the two layers close to the plates have a relative rotation of tt/2. At high temperature, in 
the isotropic phase, the boundary conditions only affect the system close to the boundaries, but when the crystal is 
quenched below the transition temperature (also called the clearing point), deep into the nematic phase, the alignment 
of the molecules with the boundary conditions on the plates extends into the bulk. In this way two different states 
appear, corresponding to the possibility of the molecules to rotate between the directions imposed by the two boundary 
plates in either a clockwise or an anti-clockwise sense. In other words, after the quench the TNLC cell develops two 
equivalent states, which we may call left-handed and right-handed. Domains of the two states are separated by 
disclination lines ||lo[| , defined as the points where the sense of rotation changes sign. The system is effectively two- 
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dimensional and the dynamics of the left and right-handed domains is very well described by nonconserved coarsening 
dynamics. 

In order to reproduce the situation studied in the present paper, it is necessary to shear the TNLC cell in a such 
a way that the shear direction is parallel to the two plates (the flow direction is of course parallel to them). Namely, 
the mutual orientation of the two plates must not be changed in the experiment, while the orthogonal walls must be 
moved in order to create the shear. In this way our (xy) plane would be parallel to the rubbed glass plates. 

Given that a vital condition for testing our asymptotic results, in the case of time-independent shear, is the 
possibility to shear the system for a long time, it seem to us that a linear geometry is probably unsuitable for such an 
experiment. On the contrary, a circular setup may be more convenient: by taking two circular glass plates, rubbed 
tangentially and radially, it is possible to create a cell whose wall, orthogonal to the plates, can now be rotated 
indefinitely. In order to create the shear it is necessary to place a fixed cylinder at the center of the system. In this 
way the material in contact with this cylinder is stationary, while the layers close to the outer walls move with a given 
tangential velocity uq, creating a velocity profile given by 

upRo f ^ Rl 

^0 



<r) - ^2 Vr2 ( - — ) ' < r < i?o , (140) 



where i?o and i?c are the radius of the cell and of the internal cylinder, respectively. If Rc — Rq Rq, it is possible to 
produce a flow identical to the one studied in the present work and to study the long time dynamics of the domains 
under shear. Indeed, by setting r — Rc + y, we have 



to be compared with relation ([l^). 

Finally, testing our results in the case of oscillatory shear should be easier from the experimental point of view, 
since the periodicity of the shear function allows for the simpler linear geometry. As we have seen, the main growth 
follows a t^/^ law, modulated by some oscillations in the longitudinal direction. In particular, it should not be difficult 
to test whether the ratio of perpendicular and parallel domain sizes satisfies the relation 



in the regime where 7 3> w. 



^ - - , (142) 

h 1 



IX. CONCLUSIONS 



In this paper we have analytically studied the effect of a shear flow on phase ordering, for a statistical system 
with nonconserved scalar order parameter. We have developed a self-consistent anisotropic version of the OJK 
approximation, by means of which we have calculated the growth exponents for time-independent shear in two and 



three dimensions (relations (|6^) and (127|), and we ha ve found the scaling form of the equal-time two-point correlation 
function in both cases (relations (||), (|7^) and (|l3|)). While for d = 3 our results are consistent with some simple 



scaling arguments and with the results obtained for conserved dynamics in the limit of large dimension N of the order 
parameter, in d = 2 we find that domain growth is so heavily affected by the shear, that the domains experience a 
narrowing which in principle makes their thickness vanish in the limit t — > 00. However, as we have pointed out, 
our calculation is likely to break down for very long times, when the interface and the domain thickness are of the 
same order. What happens beyond this stage is still unclear to us: it is possible that a time-dependent steady state 
develops, with very narrow domains coalescing and giving rise to new thicker domains, which start narrowing again. 
Another possible scenario is that when ~ ^ domains start breaking and stretching again, giving rise to a steady 
state like the one depicted in |^ . Further work is needed to clarify this point and it is to be hoped that experiments on 
twisted nematic liquid crystals, as described in the last Section, will lead to a deeper understanding of this problem. 

We have also studied the case of an oscillatory shear in two dimensions, finding a standard t^/^ growth, modulated 
by periodic oscillations which occur in opposition of phase for the parallel and perpendicular direction. Interestingly 
enough, all our results in this case arc largely independent on the particular form of the shear rate oscillations. 

It is important to note that, in two dimensions, our results satisfy the topological constraint on the growth of the 
scale area, both in the time-independent and oscillatory case. This fact, together with the results of our numerical 
simulations, strongly support the validity of our method in the study of coarsening systems under shear. 
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Of course, it would be very interesting to know whether some of our results (in particular in dimension two) are 
preserved for conserved dynamics, which is the relevant case for describing spinodal decomposition in binary fluids. 
Unfortunately, the OJK approximation cannot be used in this case, since the very starting point, the AUen-Cahn 
equation for the interface motion, does not hold when the order parameter is conserved. It is therefore still unclear 
how to go beyond the large- A'^ limit in the context of spinodal decomposition under shear. 
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